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We show that, in the continuum limit, the dynamics of Hamiltonian systems defined on a lattice 
with long-range couplings is well described by the Vlasov equation. This equation can be linearized 
around the homogeneous state and a dispersion relation, that depends explicitly on the Fourier 
modes of the lattice, can be derived. This allows one to compute the stability thresholds of the 
homogeneous state, which turn out to depend on the mode number. When this state is unstable, 
the growth rates are also function of the mode number. Explicit calculations are performed for the 
a-HMF model with < a < 1, for which the mean-field mode is always found to dominate the ex- 
ponential growth. The theoretical predictions are successfully compared with numerical simulations 
performed on a finite lattice. 



1. INTRODUCTION 

Long-range interactions are found in many different 
domains of physics: self-gravitating systems, unscreened 
Coulomb systems, plasmas, cold atoms, two-dimensional 
fluids [H-Hl- Systems with long-range interactions have 
interesting properties at equilibrium, such as ensemble 
inequivalence [6|. In this paper we will be concerned 
with their dynamical out-of-equilibrium behavior. 

In the continuum limit, where the number of parti- 
cles, N, goes to infinity, the dynamics of these systems 
is well described by the Vlasov equation @, Q- It nas 
been found that stable stationary states of this equa- 
tion correspond to quasi-stationary states (QSS) of the 
finite- N system Q . Long-range systems generically show 
a rapid evolution towards QSS, where the system re- 
mains trapped for a time that increases algebraically with 
N. In the long time limit the system eventually relaxes 
to Boltzmann-Gibbs equilibrium. This complex relax- 
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ation process was first described in pioneering papers of 
Hcnon and Lynden-Bell [ll|. 

The Vlasov equation is usually derived for Hamiltoni- 
ans of the form 

N 2 N 
3=1 jjik=l 

where denotes the position of the j-th particle in three- 
dimensional space, pj its conjugate momentum and V(x) 
the potential function. 

In this paper, we will show that systems defined on a 
one-dimensional lattice with N sites can also be described 
by a Vlasov equation, if the coupling constant decays 
algebraically with the distance between two sites of the 
lattice. 

We consider Hamiltonians of the following form 

N r,2 1 N I \ 

^2 + 2N .f- Jsi-Sfcf ( ) 

j = l 3,k=\ J 

where Xj — ja is the fixed coordinate of the j-th site on a 
one-dimensional lattice (the lattice constant a will be set 

to unity) and N = Y,j=ij~ a ~ ^ 1_Q , with < a < 1, 
is a normalization constant which makes the potential 
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energy extensive in TV [17| . We choose periodic boundary 
conditions, x x = £jv+i,qi = qw+i, Pi = Pw+i- 

The vector variable q., has now a different physical in- 
terpretation: it represents the internal degrees of freedom 
attached to site j. These degrees of freedom are assumed 
to be still described by a Hamiltonian structure with a 
momentum vector pj conjugate to q^. The dimension D 
of these vectors is not the dimension of the physical space 
in which the motion of the particles takes place, but the 
dimension of the space of internal degrees freedom, so 
that it can take values different from D = 3. 

Models of this kind were introduced by Dyson [12, [H[ , 
Thouless [lH in the 60's with the aim of understand- 
ing phase transitions in one dimension, and were also 
connected to the physics of Kondo effect [H, EH ■ How- 
ever, the internal degrees of freedom were discrete Ising 
variables while, following [l8j], we here generalize the 
model to continuous variables representing, e.g., XY 
spins. Moreover, the addition of a kinetic term in the 
Hamiltonian allows us to treat the microcanonical dy- 
namics. 

To be concrete, let us consider the case in which we 
attach to site j a two-dimensional vector pointing to the 
unit circle: q^ = (cos qj , sin qj ) , with — 7r < q,j < n an 
angle. This vector physically represents an XY rotor 
with conjugate angular momentum pj pointing towards 
the indirection. For example, let the function v(x,y) 
be the scalar product of the two arguments f(qj,qfc) = 
— Jqj q k = — J cos(qj — <?&), where the coupling constant 
J can be positive (ferromagnetic coupling) or negative 
(antifcrromagnetic coupling). This coupling has been al- 
ready considered in the literature [18l42l| and the corre- 
sponding model goes under the name of a-HMF, where 
the acronym stands for Hamiltonian Mean Field. Indeed, 
when a — the model reduces to the HMF model p^ - 
|24| : the lattice structure is removed in this limit and the 
rotors are coupled with equal strength. 

Two kinds of results were obtained for the ferromag- 
netic a-HMF model. On the one hand, it was shown an- 
alytically that the statistical equilibrium is the same as 
the one of the all-equally-coupled a = case, both in the 
canonical and in the microcanonical ensemble (2(J l2lj . 
These results were anticipated by numerical simulations 
performed in the microcanical ensemble [ID, [l9[ . On the 
other hand, simulations performed in out-of-equilibrium 
conditions have shown that no significant difference is 
observed in the features of the QSS with respect to the 
a = case [TfJ [25|, [26| . This raises the question of the 
full equivalence between the model with all-to-all cou- 
pling and the lattice model (see also (27|). 

We will here show that the Vlasov formalism can be 
actually generalized to long-range interactions on a lat- 
tice. In the TV — > oo limit, each infinitesimal element of 
the lattice contains an infinity of sites, to which a local 
distribution function in phase-space /(q, p; x, t), depend- 
ing also on the spatial coordinate x along the lattice, can 
be associated. A Vlasov equation for this local distri- 
bution function can be then derived, and the long-range 



interaction determines a global coupling between these 
local distributions at different spatial locations. A lin- 
ear stability analysis of a state which is both uniform 
in the internal degrees of freedom and along the lattice 
will be performed and the stability thresholds will be de- 
termined analytically. When this state is unstable, dif- 
ferent Fourier modes grow with different rates and the 
zero mode, i.e. the mean-field mode of the a = model, 
is found to dominate the instability process, indicating 
that, indeed, the dynamics shows similar features with 
the HMF model. 

This paper is organized as follows: Section[5]is devoted 
to the derivation of the Vlasov equation for long-range 
interacting systems on a lattice. In Section [3] we discuss 
the linearization around the uniform state and we find a 
dispersion relation expressed in terms of lattice Fourier 
modes. In Section|4]wc derive the stability thresholds and 
the growth rates of the Fourier modes for the a-HMF 
model, which are tested against numerical simulations 
performed with finite TV in Section[5] Finally, in Section[5] 
we draw some conclusions. 



2. VLASOV EQUATION ON A LATTICE 

In order to simplify the derivation, we restrict to one- 
dimensional internal degrees of freedom attached to the 
lattice site j . These degrees of freedom are represented by 
the canonically conjugate pair (qj,pj). The Hamiltonian 
is the following 



H 
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■x k \ a 



(3) 



where v(x,y) is symmetric in the two arguments. The 
equations of motion are then 



9j 
Pj 



Pj- 
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\Xj-X k \ a 



(4) 
(5) 



where the prime denotes the derivative with respect to 
one of the two arguments. Since we have chosen periodic 
boundary conditions, we use the closest distance conven- 
tion, i.e. the distance \xj — x k \ actually corresponds to 
mm(\xj - x k \, 1 - {\xj - x k \)). 

In the continuum limit, the lattice (which for periodic 
boundary conditions can be thought to be a circle) is 
densely filled with sites: each infinitesimal lattice ele- 
ment dx contains a diverging number of sites as TV goes 
to infinity. Therefore, one can define in this limit a distri- 
bution function f(q,p;x,t), which depends both on the 
internal variables {q,p) and on the coordinate x along 
the lattice, which is confined to the interval [0, 1]. In the 
continuum limit, the equations of motion become 



<T 
P 



Pi 



(6) 

f( q ',p'. x ',t)0^dq'dp'dx', (7) 
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where IS such that 1/JV 



JV-s-oo 



> K a /N. The latter 

-1/2 



limit can be shown to yield k 1 = f_ 1 j~ dx/\x\ a . Equa- 
tion Q naturally leads to the definition of the long-range 
potential 

V x [f]{q,t) = n a fffd q 'dp'dx'f(q',p';x',t)- V( 



\x - x'\ c 



(8) 

which determines the motion of the conjugate variables 
(<7,p) at point x. This definition allows us to write the 
following Vlasov equation 



(9) 



where the prime refers to d/dq. It is interesting to note 
that this equation also derives from the following Hamil- 
tonian 



H[f] = 



(10) 

Thus, Eq. ([9]) actually represents a continuum of Vlasov 
equations along the lattice, globally coupled through the 
potential (JSJ. The absence of a spatial derivative d x in 
Eq. ((9]) reminds us that the lattice sites are fixed. 

More rigorously, one should follow a procedure that 
begins with the definition of the so-called empirical mea- 



N ^ 
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f e (q,p;x,t) = — y^6(q - qj(t))S(p - pj(t))S(x - Xj) , 

(11) 

and then derive the Klimontovich equation as done in 
[35l [36j . The Vlasov equation could be then derived using 
a convenient averaging method. 

It is also important to stress that definition ((5J is con- 
sistent only for a < 1, since the integral J dx'/\x — x'\ a is 
convergent around x only for such values of a. For a = 1 
the integral diverges logarithmically, we will not consider 
this case. 

Let us comment about the normalization of /, an is- 
sue linked to the density of sites along the lattice. Since 
we are here considering a homogeneous lattice, each in- 
finitesimal element dx contains the same number of par- 
ticles. Hence the norm of f(q,p;x,t) has to be the same 
at any position x 



f(q,p;x,t)dqdp = 1. 



(12) 



Yet, our formalism is also suitable for heterogeneous lat- 
tices with non- uniform density of sites p(x), for which 
the normalization is given by 



f(q,p;x,t)dqdp = p(x). 



(13) 



The density p(x) enters as a weight when deriving the 
local potential V^[/]. 



Finally, we note that the generalization to higher- 
dimensional lattices and to higher dimensional spaces of 
internal degrees of freedom is straightforward. For D- 
dimensional conjugate variables (q,p) defined at posi- 
tion r on a d-dimcnsional lattice, the l/\xj — Xk\ a cou- 
pling turns into l/\rj— Tfe| a , and the distribution function 
f(q,p;r,t) obeys a Vlasov equation at each site r. 



3. THE DISPERSION RELATION 

As anticipated in the introduction, the emergence of 
a Vlasov equation in long-range systems is of particular 
importance for their dynamics, since it is responsible for 
the long-lasting out-of-cquilibrium regimes for finite N, 
the so-called QSS. This phenomenon is however bounded 
to the existence of Vlasov stationary stable states. Wc 
shall now turn to determine the conditions for such states 
to be present for the Vlasov equation on a lattice © . Let 
us consider a stationary solution fo of Eq. §§§ 



pd q fo-V^[f ](q)d p f = 0. 



(14) 



From the above relation, it appears that the motion of 
a test-particle at position x can be derived from the fol- 
lowing Hamiltonian 



h x (q,p) 



EL 
2 



V x [fo](q) 



(15) 



Following the observation that h x (q,p) is a conserved 
quantity for stationary solutions /o, several authors [28l — 
HH H3 studied inhomogeneous Vlasov equilibria express- 
ing the stability problem in action-angle variables. How- 
ever, for the sake of simplicity, wc will here focus on 
homogeneous stationary states, both in position x along 
the lattice and in the internal variable q. This implies 
that f = f (p) and V^[f ](q) = for all x. The stability 
of /o can be studied by considering a small perturbation 
of the stationary distribution 



f{q,p;x,t) = f (p) +6f(q,p;x,t). 



(16) 



Inserting the above expression into the Vlasov equation 
and dropping the quadratic term in Sf leads to the 
linearized Vlasov equation 

d t (Sf)+pd q (6f) - MV'lSfUt) = 0. (17) 

If we now focus on an cigenmode 8ft{q,p]x) = 



f(q,p; x) of the linearized dynamics, we get 



(X+pd g )f(q >P ;x) - MV^[f}(q) 
which can be rewritten as 

A — 



0. 



(18) 



e x lf 



P 



-fo(p)v:if](q) = o- (i9) 



The integration over q yields the dispersion relation 



f(q,F, x) - fo(p)- 



e xq *V:[f](q'W = 0, (20) 
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where we have assumed the integration constant to be 
zero. The solution of this dispersion relation depends on 
the form of the local potential V x [/] . If we restrict to the 
derivation of the stability threshold, we can set A = in 
the previous expression and perform the integral over q 1 ', 
which gives 



V x [f](q,p;x) = 0. (21) 



This is indeed a set of equations for each position x, all 
coupled together by the V x [f] term. 

Now, since we have considered a lattice with peri- 
odic boundary conditions, it is natural to introduce the 
Fourier operator 

&k[f](q,p) = ^J dx e- 2iwkx f(q,p;x) = f k (q,p), 

(22) 

and rewrite the distribution / as a sum of Fourier modes 
f(q,p;x) = ^f k (q,p)e 2i * k *. (23) 



The potential V x [f] exhibits the interesting property of 
being diagonal in Fourier space [391 ] . Indeed, inserting 
the sum (12"31 into the definition JHI results into 



V x [f](q) = n a J2 ( / dx ' 



\x - x'\ c 



(24) 



X // dq'dp' f k (q',p')v(q,q'). 



Recalling the definition of the nearest distance on the pe- 
riodic lattice, the above integral over x' is actually per- 
formed from x — 1/2 to a; + 1/2, so that it is appropriate 
to introduce the change of variables x' — > y = x' — x. 
One gets 



x+1/2 2iirkx' 

dx'- 



x-1/2 \ x ' — X Y 



where c k {a) is defined as 



/+1/2 „2inky 
-i/2 \y\ a 



Jli-nkx 



f +1 / 2 dy 



Ck{a) = k c 



+ 1/2 e 2i-Kky 



1/2 \y\ a 



dy = 



+ 1/2 e 2inky 



1/2 \y\ 



dy 



-1/2 



(26) 



1/2 \y\ a 



dy 



Equation (|24|) leads us to define the mean-field poten- 
tial V(q) for a distribution g(q,p), which does not depend 
neither on a nor on the position on the lattice 

V[g](q) =11 dq'dp'g(q',p')v(q 7 q') . (27) 



Remark that in the limit a = 0, when the long-range 
potential V x in formula ([8]) does not depend any more 
on the position x along the lattice, one has V[f] = 
jj dq'dp 1 f{q' ,p')v(q, <?'), exactly as in mean-field mod- 
els. Equation ([27)) allows us to rewrite the dispersion 
relation (f2"U|) as 



-\3. rq 



:/o(p) 



P 



e*TV'[f k ](q')d q '.(28) 



Thus, applying the inverse Fourier transform leads to the 
following set of decoupled dispersion relations for each 
Fourier mode f k 

f k (q,p) - Ck (a)f (p)'—^ f e X ^V'[f k ](q')\dq'\ = 0. 



(29) 

The quantity Afc represents the eigenvalue associated to 
mode k. It is interesting to note that, although the so- 
lution of these equations will generically depend on the 
specific form of the potential V[f], the dependence on 
the specific Fourier mode only appears through the coef- 
ficient c k (a). We shall see in the next Section that this 
allows us to derive the growth rate for all Fourier modes 
in the case of the a-HMF model. 



4. STABILITY THRESHOLDS AND GROWTH 
RATES FOR THE a-HMF MODEL 

As a test-bed for systems of rotators on a ring-lattice 
with long-range couplings, we now turn to analyzing the 
dispersion relations for each Fourier mode of the a-HMF 
model [l8l42~ll ]. The rotators are thus coupled through 
a two-body potential v(q,q') = — cos(q — q'), which is 
tuned down by a decay over the distance along the lattice 
l/\x - x'\ a . The potential V, defined in Eq. ([27]) reads 



V[f](q) = -Kc, [M x [f] cosq + M y [f] sinq) , (30) 

where the magnetization M[f] is defined as 
M[/] = M x [f]+iM y [f] (31) 
dq'dp' f(q',p') cosq' + i J j dq'dp' f{q',p')smq'. 

Equation (f29j) turns into 



= fk(q,p) ~ c k (a)M 



(32) 



A — 

x / e k p 



M x [f k ]sinq' -M y [f k ] cos q') \dq' 



This equation can be solved by multiplying both sides 
by either cos q or sin q and integrating over the phase- 
space (<?, p) , which leads to a linear system of equations 
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in M x [f k ] and M y [f k ] 



Thus, the eigenvalue of the fc-th Fourier mode is given by 



M x [fk] 1 - c fc (a)4VI/o] + M y [f}c k (a)I^ x [f a ] = 



M x [f k ]c k (a)I^[f ] - M y [f] (1 + c k (a)I^ x [f a ] ) = 0(33) 

where I x Y [fo] is defined, for a scalar A and two functions 
X(q),Y(q) (here X(q) = cosq and Y(q) = sing), as 



dp^M I & qe - X ix{q) [ dq'e^Ylq'). 
P J J qo 

(34) 

Assuming that cigcnmodes / have a non-zero magnetiza- 
tion, the system of equations (|33j) implies 



= (l - c k (a)I% Y [fo]) (l + c k (a)I^ x [f 

+c k (a) 2 I^[fo}I^ x [fo}- (35) 



Using the following formulae 

f'l 



I „A fc 



dq' e 



p smq = 



dq' e Xk p cos q' = 



, 2 f — sin q - cos q I (36) 
1 + V P 

p 2 



1 



At 

p 2 



smq 



V 



■ cos q I (37) 



one can calculate the integrals ([34]). so that dispersion 
relation ([35"]) can be finally rewritten as 



= \l + nc k (a) [ dp^ 

v J K 1 



+ Trc k (a) / dp 




(38) 



Note that in this relation only the coefficient c k (a) 
changes from one Fourier mode to another. 
Let us now consider a waterbag distribution 0| 



hip) 



i i 



(Q(p + Ap)-Q(p-A P )) 



(39) 



2tt 2Ap 

where is the Heavyside function. It's first derivative is 
1 1 



2tt 2Ap 



(S(p + Ap)-S(p-Ap)) , (40) 



where S is the Dirac delta function. Since this distribu- 
tion is symmetric in p, the last term in Eq. (|38[) vanishes, 
and the dispersion relation simplifies into 



= 1+ TTCfe 



(a) ( dp—J- 
J v 1 



Cfc(a) 



(41) 



c k {a) 



Ap 2 



(42) 



Hence, mode k is stable provided Ap 2 > c k (a)/2. If 
instead Ap 2 < c k (a)/2 the mode is unstable and grows 
with the rate Rc(Afe). 



5. NUMERICAL SIMULATIONS 

We shall now compare the predictions of the previous 
Section concerning the stability of the waterbag initial 
state for the a-HMF model with numerical simulations. 
We integrate the equations of motion (|4I5[) using an opti- 
mized fourth-order symplectic scheme [38| with time step 
0.1, which guarantees a good conservation of the energy 
during time evolution. We do not need here to perform 
long-time evolutions, because we have just to look at the 
initial growth of the Fourier modes. Hence, we can push 
as much as we can towards larger values of N: we report 
here simulations with N up to 2 18 . The reason for choos- 
ing a power of 2 for N is that we use a Fast Fourier Trans- 
form algorithm in order to compute the force in Fourier 
space: this allows us to perform simulations in O(NmN) 
time instead of the 0(N 2 ) time needed to compute the 
double sum over the sites of the lattice. We recall that 
for the a-HMF model v(qj,q k ) = — cos(qj — q k ) and we 
choose here a = 0.8. 

The system is initiated by spreading the particles ran- 
domly in the rectangle [0; 2ir] x [—Ap; Ap] in the phase 
space (q,p) attached to each site of the lattice. The time 
evolution of each harmonic of the magnetization is fol- 
lowed using the observable 



m k 



—ikx i 

e e 



f(q,p;x,t)dq dp 



1 

N 



E 



(43) 



(44) 



The simulations allow us to observe an exponential 
growth of m k (t) (see Fig. QJa,b)), with a rate that de- 
pends on the harmonic number k, as expected from for- 
mula (|4"2"j) . They also confirm that some higher harmonics 
may be stable while lower ones are unstable: for example, 
for Ap = 0.65, formula (|42j) predicts that Fourier modes 
with k > 1 should be stable. Simulations indeed reveal 
that for such a value of Ap, mode k = 3 does not grow 
exponentially any more (see Fig. HJb)). 

Let us give some further informations about the time 
evolution of m k (t). The harmonics are initially at a level 
of order 1 / v^V, as expected from a random distribution. 
The very short-time dynamics (t < 2 for Ap = 0) does 
not allow us to observe the exponential growth, presum- 
ably because there is still a strong competition between 
the numerous eigenmodes, be they stable or not. The 
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most unstable mode eventually dominates, and it is ob- 
served that the exponential growth sets in for the time 
range 2 < t < 7 for Ap = 0, which is almost the same for 
all the harmonics (the growth rate is calculated over this 
time window). Remark that for Ap = 0.65, in Fig. [TJb), 
although the growth is very close to an exponential one 
for mode k = 0, it is not for k = 1,2,3, as predicted by 
the theory. 

Finally, the results of a simulation with much less par- 
ticles TV = 2 12 are drawn in Fig. Qlc), to illustrate the 
finite-size effects. In this case the Fourier modes start for 
t = at a much higher level, so that the growth occurs 
on a much shorter time window. It becomes then very 
hard to identify the exponential character of the growth. 

In order to test more thoroughly the accuracy of the 
theoretical predictions, simulations for homogeneous wa- 
terbags with different values of Ap were performed. In 
Fig. [2] we show the growth rates as a function of Ap for 
different Fourier modes. The agreement between the- 
ory and simulations is excellent for the mean-field mode 
(k = 0), quite good for the two following modes k = 1, 2 
but it gets worse as the harmonic number increases, es- 
pecially for the larger values of Ap. This phenomenon 
is possibly due to the fact that higher Fourier modes are 
more sensitive to finite-size effects. Finally, let us com- 
ment on the choice a = 0.8 for our test: as can be eval- 
uated from Eq. p^j) . the domains of instability of k > 1 
modes is very narrow for smaller values of a, so that a 
a close to unity was best to present a stability diagram 
showing clearly the successive extinction of the Fourier 
modes. 



6. CONCLUSIONS 

The main result of this paper is the derivation of the 
Vlasov equation for a lattice with long-range couplings. 
As an application, we have studied the stability of the 
linearized Vlasov equation for the a-HMF model, deriv- 
ing a dispersion relation which is expressed in terms of 
Fourier modes . 

We have computed the stability thresholds for the dif- 
ferent Fourier modes and the growth rates when the 
modes are unstable. We have shown that the fastest 
growth rate is the one corresponding to the mean-field 
zero mode. This agrees with the widely spread opinion 
that the a-HMF model with < a < 1 behaves in all 
respects as the HMF model, which corresponds to the 
limit a = 0. 

In order to support this conjecture, we have also tried 
to identify states where higher harmonics would domi- 
nate, performing numerical simulations in which the sys- 
tem is initialized in a state which is spatially modulated 
with a given wave-number in order to promote the time 
evolution of higher modes. In all cases, we have observed 
the growth of other modes, including the zero mode. This 
latter eventually takes the lead and dominates over the 
others. 
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Figure 1: Time evolution of the Fourier modes of the magneti- 
zation mk(t) (|44[) for the a-HMF model. The initial waterbag 
distribution is uniform in q £ [0, 2n] and in p 6 [—Ap, Ap] 
with Ap = in panel (a) and Ap = 0.65 in panels (b) and 
(c). For these values of Ap the homogeneous waterbag is un- 
stable. The Fourier modes rrih are plotted vs. time for k = 
(thick full line), k = 1 (thick dash-dotted line), k = 2 (thick 
dashed line) and k = 3 (thick dotted line). The thin lines are 
exponential fits. The two upper panels have N = 2 18 , while 
the lower panel has N = 2 12 and shows the strong finite-size 
effects: the exponential growth is no more clearly visible. 



A more specific remarks concerns the amplitude Cfc(a 
of the Fourier modes of the potential, given in formula 
(|2l)l) . In the limit a = all the harmonics k > 1 van- 
ish, giving the potential in the mean-field limit, which 
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Figure 2: Growth rates Re(Xk) of the first Fourier modes of 
the magnetization as a function of Ap. The initial conditions 
are homogeneous waterbags. The lines correspond to the the- 
oretical prediction (|42|l for k = (full),fc = 1 (dash-dotted), 
k = 2 (dashed) and fc = 3 (dotted). The results of numerical 
simulations are marked by the symbols (resp. o, *, x and □). 
For each value of Ap ten runs were performed: the symbols 
represent the average, while the error bars stand for the min- 
imal and maximal values. The symbols for k = 1 and k = 3 
are slightly shifted in Ap for the sake of clarity. 



corresponds to the one of the HMF model. On the other 
hand, when going toward the limit a — > 1 from below, one 
gets Cfc(a) = 1 for all fc's. This means that all the Fourier 
modes, and thus all the length scales, become equally rel- 
evant in this limit. For a = 1, the series of Fourier modes 
that represents the potential is no more summable and, 
as a consequence, the potential (J5J) is ill-defined: it's the 
signature that one is leaving the long-range region. 
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